Transverse momentum dependence of Jjip shadowing effects 
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We present a new approach to estimate the effect of the gluon shadowing in nucleus+nucleus 
collisions and its consequences on the J/tp production yield. Using kinematical information available 
from the measured J/tp production in proton+proton collisions at ^/sNN ~ 200 GeV, we build a 
Glauber Monte Carlo code which takes into account shadowing in two alternative ways: multiple 
scattering corrections or evolution of parton densities. We exploit the dependence of these 
different parameterizations to the J/ip transverse momentum and we give the first predictions on 
the resulting pr dependence of the nuclear modification factor in deuteron+gold collisions at the 
same energy. 

PACS numbers: 25.75.Dw 25.75.-q 24.85.+p 



The transition from hadronic matter to a deconfined 
state of quarks and gluons, the so-called Quark-Gluon 
Plasma (QGP), is the subject of intense experimental and 
theoretical works in heavy ion physics. Recent results on 
J/^ production measured by the PHENIX experiment at 
the BNL Relativistic Heavy Ion Collider show a signif- 
icant suppression of the J/^ yield in Au+Au collisions 
at ^/sNN— 200 GeV [H- A possible origin of this sup- 
pression is the screening of the charmonium pair in the 
hot and dense nuclear medium, as expected in case of a 
QGP formation. Nevertheless, PHENIX data on d-|-Au 
collisions have also revealed that cold nuclear matter 
effects play an essential role at these energies. In fact, 
in order to determine the importance of hot and dense 
matter effects on J/tp production in ^-t-A collisions, it is 
fundamental to have a good understanding and a proper 
baseline for cold nuclear matter effects including initial- 
state shadowing corrections. 

Motivated by this goal, wc have developed a Monte 
Carlo code, based on the Glauber model, where the 
cold nuclear matter effects - shadowing - have been 
introduced. The inputs of our Monte Carlo are the 
J/V' rapidity and transverse momentum distributions ex- 
tracted from y/sNN— 200 GeV p + p data measured by 
the PHENIX experiment We investigate two differ- 
ent possibilities for the shadowing corrections: on one 
side we consider a multiple scattering approach, where 
the coherence effects lead to a modification of the nu- 
clear ratio; on the other side we apply a model based on 
Q^-evolution of the nuclear ratios of parton distribution 
functions (PDF). 

We present our results for both models on the nuclear 
modification factor RdAu as a function of rapidity y, num- 
ber of collisions A^coU and transverse momentum pT- For 
the first time, we present the px dependence of the shad- 
owing effects on the J/ip- We found significant differences 
according to the type of shadowing considered, hence giv- 
ing a new opportunity to test shadowing schemes. Note 
that, in order to show up shadowing effects, we have not 



included the contribution from the nuclear absorption, 
i.e. the J/tp break-up by nucleons in the final state. 

Shadowing refers to the mechanism that makes the nu- 
clear structure functions in nuclei different from the su- 
perposition of those of their constituents nucleons. Sev- 
eral explanations have been proposed. Here, we will con- 
sider two different approaches : 

• On one hand, in the rest frame of the nucleus, 
nuclear shadowing can be understood as a con- 
sequence of multiple scattering 0, the incom- 
ing virtual photon splitting into a colorless qq 
dipole which interacts with the nucleus with typ- 
ical hadronic cross sections resulting in absorption. 

• On the other hand we will consider the evolu- 
tion models These approaches provide the 
PDF at a given initial value of and study its 
evolution through the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) equations. 

Following the first approach - models based on multiple 
scattering - Q , one can calculate the reduction factor due 
to the shadowing corrections at fixed impact parameter 
b for a nucleus as 
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The fmiction F{y,pT) is given by the integral of the ra- 
tio of the triple Pomeron cross-section over the single 
Pomeron one - equivalent to the reduction due to the 
interaction among the gluons - 
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where ymin = In ( ^•^'" ) and ymax = | In TV- y 

is the center of mass rapidity of the produced particle, 
y > for the projectile hemisphere and y < for the 
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target one. tut is the transverse mass, rriT = \/ 
The parameters C = 0.31 fm^ and A = 0.13 are fixed 
through data on DIS scattering In the following we 
will refer to this model as CF shadowing. 

The second type of models study the Q^-evolution of 
nuclear ratios of PDF, 



Rt{x,Q^) = 



nucieon 



(x,Q2) ' 



fi = q,q,9, (3) 



through the DGLAP evolution equations. Nuclear ratios 
are parameterized at some initial scale Qq ^ 1 ~2 GeV^ 
which is assumed large enough for perturbative DGLAP 
evolution to be applied reliably. These initial parame- 
terizations for every PDF have to cover the full x range: 
< X < 1. Then these initial conditions are evolved 
through the DGLAP equations towards larger values of 
and compared with experimental data. From this 
comparison the initial parameterizations are adjusted. 
The centrality dependence of shadowing is not addressed 
in these models as the existing experimental data do 
not allow its precise determination, although some ap- 
proaches provide an ansatz for such a dependence. It can 
be parameterized [7] assuming that the inhomogeneous 
shadowing is proportional to the local density, z), 



Rf{b,x,Q') = l + [Rf{x,Q'' 



JdzpA{b,z) 
J dzpA{0,z) 



where b and z are the transverse and longitudinal lo- 
cation in position space, pAib,z) corresponds to the 
Woods-Saxon distribution for the nucieon density in the 
nucleus, related to the nuclear profile function Ta(6) 
by / dzpA{b,z) = A TA{b), po is the central density, 
given by the normalization / d^b J dzpA{b,z) = A and 
Rf{x, Q"^) is the shadowing function from EKS98 as de- 
fined by eq. The integral over z includes the mate- 
rial traversed by the incident nucieon, so we are con- 
sidering that the incident parton interacts coherently 
with all the target partons along its path length. Np 
is the normalization factor, and it is chosen so that 
{l/A) ^ d'^b ^ dzpA{b,z)Rt{b,x,Q'') ^ Rt{x,Q''). In the 
following we will refer to this model as EKS shadowing. 

Within our framework, we study the modification of 
the J /ip production in nuclear matter. It is often pre- 
sented using the nuclear modification factor: 



Rab — 



dN 



AB 



{N,ou)dN^l^ 



(5) 



where dN^/g and dNpl^ are respectively the J/-0 yield 
observed in ^ -I- i? collisions and va. p + p collisions, and 
(Ncoii) is the average number of nucleon-nucleon colli- 
sions occurring in a, A + B collision. For a hard process, 
such as J/tp production, in absence of nuclear effects, the 
Rab ratio should equal unity. 



In the present study, this quantity is obtained with 
a Glauber Monte Carlo. Within this code we describe 
the collision of two nuclei, determine the number of nu- 
cleons in the path of each incoming nucieon and calcu- 
late the number of nucleon-nucleon collisions. Nuclear 
density profiles have been defined using a Wood-Saxon 
parametrization for any nucleus A > 2 and the Hulthcn 
wavefunction for the deuteron The nucleon-nucleon 
inelastic cross-section at y^sjvw — 200 GeV has been 
taken to 42 mb and the average nucieon density to 0.17 
nucleons/fm'^. 

In absence of nuclear effects, the J/tp yield in a nucleus- 
nucleus collision corresponds to a simple superposition 
of the yield obtained in p + p collisions multiplied by 
the average number of nucleon-nucleon collisions. In 
p + p collisions, the J/ip production has been measured 
at y/sNN = 200 GeV by the PHENIX experiment, pro- 
viding the J/ip differential cross section versus both the 
transverse momentum px and the rapidity y Q . 

Based on the parameterizations used to fit these two 
distributions we can generate a sample of J/tp'- we ran- 
domly pick y, Pt and (p (the azimuthal angle, described 
as a flat distribution within [0,27r]) and use these val- 
ues as the inputs of our Monte Carlo. Then, using the 
shadowing parametrization, we can compute R^^^^w ' 
nuclear modification factor when the J/tJj production is 
affected by shadowing. 

Figure[T]shows the distributions used in our code. Note 
that in order to restrain ourself to a physical phase space 
domain, we require for each J/ip the kinematic conditions 
< xi,X2 < 1 where 



rriT 
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\/snn 



and X2 = 



y/SNN 



(6) 





Pj (GeV/c) 



FIG. 1: J/ip rapidity and pr distributions in p + p collisions 
extracted from Q]. Vertical scales are arbitrary units. 



As previously mentioned, the shadowing models pre- 
dict that the modification of the structure functions de- 
pends on the nuclear environment surrounding the in- 
teracting partons. They both parameterize this spatial 
dependence as a function of A Ta(6) where A is the nu- 
cleus mass number and Ta(6) is the nuclear profile func- 
tion at fixed impact parameter &, normalized to unity. 
In fact, A TA{b) corresponds to the number of nucle- 
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ons seen by a given nucleon at a distance b. This quan- 
tity depends on the position of the given nucleon in the 
nucleus; central nucleons see more surrounding nucleons 
than peripheral ones. Within our framework, this num- 
ber can be computed event by event, thus taking into 
account the spatial nucleon distribution for the current 
nucleus. Defining as Ua the nucleon which suffers shad- 
owing and d the distance between and a nucleon Ub 
belonging to the same nucleus, tIq suffers shadowing from 
Ub if TTc?^ < atr, where atr corresponds to the nucleon 
transverse area. Since the average nucleon density in a 
nucleus is po — {^nr^)^^ =0.17 nucleons/fm^, we take 
atr = Trr^ = 3.94 fm^. We can then compute Nf^, that 
is, the number of nucleons which add a shadowing con- 
tribution to Ua. Finally, R^^^^^ is computed following 
respectively [5| and [H] for CF and EKS shadowing. 

In the context of the CF model, -R^f^ow computed 
using the Schwimmer expression given in Eq. ([1]). For 
A + B collisions, we get: 



AB 
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''shadow 



= t X - (7) 

l + NtiF{y,pT) \ + N^.F{y,pT) 

where F{y,px) is the CF shadowing function defined in 
Eq. ^ , and y and pt are the rapidity and the transverse 
momentum of the J/ip candidate. 

In the EKS calculation , the authors provide the ra- 
tio Rf of the PDF in a proton of a nucleus A to the PDF 
in the free proton. They give a numerical parametriza- 
tion of Rf{x,Q^) for all parton flavours. Here, we re- 
strain our study to gluons since, at high energy, J/ip is 
essentially produced through gluon fusion The nu- 
merical parametrization is valid for: A > 2 (otherwise 
return 1) where A is the atomic number, 10~^ < a; < 1 
and > 2.25 GeV^. Within our framework, we used 
the table provided in where the inputs are A, x and 
and the output is Rf{x, Q"^) which gives homogeneous 
shadowing. In order to include inhomogeneous shadow- 
ing, meaning that shadowing should depend on the inter- 
acting parton spatial position, we follow the prescription 



;on sp. 

introduced in [7i| and define R^^^w 
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AB 

shadow 
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-^shadow -^shadow 
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[R^{x^,Q^)-l] 
[i?f(x2,Q')-l] 



(8) 



where N^^^^/{Nt 



^^^^) reflects the local nuclear density. 
Note that we implicitly assume here that the incident 
parton interacts coherently with all the target partons 
in its path. 

The xi and X2 values are obtained from y and px as in 
Eq. Doing so is equivalent to consider the hard pro- 
cess g+g cc^ where the cc nonzero is inherited from 



the initial gluon transverse momenta. For the square of 
the momentum transfer we follow the prescription 
from 01 for the mass term and we add a contribution for 
the transverse momentum. Thus, we take: 

g2 = (2me)2 + {pT? (9) 

where ruc = 1.2 GeV/c^ is the c quark mass. 

We study the behavior of the nuclear modiflcation fac- 
tor RdAu in d-f Au collisions at ^snn — 200 GeV. RdXn is 
presented as a function of rapidity y, number of collisions 
A'coii and transverse momentum px- In both shadowing 
cases, we consider two possibilities for the J/V' trans- 
verse momentum distributions: no transverse momentum 
(pt=0) to rely on previous studies and px distribution 
extracted from the measured p+p transverse momentum 
distribution in the central {pt from |y| < 0.35) and the 
forward {pt from 1.2 < |y| < 2.2) rapidity regions. 
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FIG. 2: Nuclear modification factor as a function of rapidity 
from CF (top) and EKS (bottom) models. Three input p + p 
Pt distributions are used: pt=0, pr from |y| < 0.35 and pr 
from 1.2 < \y\ < 2.2. 

Figure [2] shows the nuclear modiflcation factor RdKu as 
a function of y for CF and EKS shadowing. In both cases 
the three px configurations are considered. A clear dif- 
ference is observed between the two shadowing models at 
the transition between the antishadowing and the shad- 
owing regions. The antishadowing region is at smaller y 
for CF shadowing (around —3 < y < —1.5) than for EKS 
shadowing (around —3 < y < —0.5). As a consequence, 
the amount of antishadowing is larger for EKS compared 
to CF. Note that CF was originally formulated to study 
the amount of suppression for the nuclear ratio due to 
coherence in the shadowing region, so its kinematic rel- 
evant limits for the present energies are fixed at around 
y > —2. The transverse momentum does not affect much 
the resulting i?£iAu since the average J/ip Pt inp + p col- 
lisions is smaller than the mass of the J/ ip (less than 2 
GeV/c). 

The nuclear modification factor RdAu plotted as a func- 
tion of the number of collisions A^coU is shown in Figure [3] 
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FIG. 3: Nuclear modification factor as a function of the num- 
ber of collisions from CF (left) and EKS (right) models. From 
up to down: backward (—2.2 < y < —1.2), central (\y\ < 0.35) 
and forward (1.2 < y < 2.2) rapidity regions. 



FIG. 4: Nuclear modification factor as a function of the J/ip 
transverse momentum from CF (left) and EKS (right) models. 
From up to down: backward (—2.2 < y < —1.2), central 
{\y\ < 0.35) and forward (1.2 < y < 2.2) rapidity regions. 



and leads to the same conclusions concerning the use of 
the real px distributions. Even though the effects are 
larger for EKS shadowing, the differences between the 
different pt configurations are pretty small and would 
probably be out of reach of experimental precision. On 
the other hand, the difference between the two shadowing 
models is clearly visible in the backward rapidity region 
(—2.2 < y < —1.2) where EKS shadowing exhibits a large 
RdAu increase while it is almost flat for the CF shadow- 
ing case. This behaviour can be directly related to the 
transition point between antishadowing and shadowing 
regions, as mentioned above. 

In contrast with previous figures. Figure [4] clearly illus- 
trates the gain in using non zero transverse momentum 
distributions instead of pt=0 as done in previous studies. 
The RdAu ratio is plotted as a function of pt. It shows a 
dependence which can reach an amplitude of about 20 %. 
The strongest pt dependence is observed for the EKS 
shadowing with a specific behavior at backward rapidity 
(—2.2 < y < —1.2) where RdAu decreases with increasing 
J/V' transverse momentum. Indeed, within EKS frame- 
work, a larger in the antishadowing region leads to a 
smaller RdAu factor. A larger pT leads to a larger Q^, 
so we expect a decrease of RdAu in the antishadowing re- 
gion and an increase otherwise. This effect is illustrated 
in Figure [2] where the use of real pT distributions induce 
a reduction of RdAu for y < —1.5 while this tendency is 
reversed for y > — 1.5. In CF shadowing, the inclusion 
of Pt in Eqs. ^ and ([2]) reduces the shadowing correc- 
tions, thus decreasing the function F{y,pT)- It results 
in an increase of the nuclear modification factor for all 
rapidity regions. 

In conclusion, we have compared, in the framework 
of a Glauber Monte Carlo, the results concerning shad- 
owing effects from two different theoretical models: CF 



and EKS. We have focused our attention on a new 
subject, the study of the shadowing dependence on 
Pt, where both models show a different qualitative be- 
haviour. In general, by comparing the two models, we 
found the amount of suppression from shadowing cor- 
rections stronger in the case of the CF multiple scatter- 
ing approach. However, when introducing the px con- 
tribution, a notable exception is found in the backward 
rapidity region where EKS shadowing show a stronger 
suppression at large pt ■ 
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